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A new data analysis method based on physical observables for WIMP dark matter searches with 
noble liquid Xe dual-phase TPCs is presented. Traditionally, the nuclear recoil energy from a scat- 
ter in the liquid target has been estimated by means of the initial prompt scintillation light (SI) 
produced at the interaction vertex. The ionization charge (C2), or its secondary scintillation (S2), 
is combined with the primary scintillation in Log(S2/Sl) vs. SI only as a discrimination parameter 
against electron recoil background. Arguments in favor of C2 as the more reliable nuclear recoil 
energy estimator than SI are presented. The new phase space of Log(Sl/C2) vs. C2 is introduced 
as more efficient for nuclear recoil acceptance and exhibiting superior energy resolution. This is 
achieved without compromising the discrimination power of the LXe TPC, nor its 3D event recon- 
struction and fiducialization capability, as is the case for analyses that exploit only the ionization 
channel. Finally, the concept of two independent energy estimators for background rejection is pre- 
sented: E2 as the primary (based on C2) and El as the secondary (based on SI). Log(El/E2) vs. 
E2 is shown to be the most appropriate phase space in which to evaluate WIMP signal candidates. 

PACS numbers: 14.80.Ly; 21.60. Ka; 29.40.Mc; 95.35.+d 
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I. INTRODUCTION 

Considerable astronomical and cosmological evidence 
supports a self-consistent ACDM model of the Universe 
where approximately 23% is in the form of dark matter. 
Discovery of the nature of this dark matter is recognized 
as one of the greatest contemporary challenges in science, 
fundamental to our understanding of the Universe. The 
most compelling candidates for dark matter are Weakly 
Interacting Massive Particles (WIMPs) that arise natu- 
rally in several models of physics beyond the Standard 
Model pQCg. 

WIMPs may be directly detected through the energy 
deposited by a recoiling nucleus following a rare WIMP 
scatter with standard matter. Given the kinematics, the 
recoil would generate an extremely small signal of <100 
keVnr (keV nuclear recoil) requiring detectors to exhibit 
good sensitivity (high detection efficiency) at low ener- 
gies (down to a few keV nr ) with good energy resolution to 
determine the WIMP energy spectrum and estimate the 
WIMP mass. Simultaneously such an ideal detector must 
have good rejection power for background events from 7- 
rays. Despite this, the experimental sensitivities of direct 
searches are now entering into the theoretically favored 
parameter space and a positive detection may be immi- 
nent from current or next generation devices. Two-phase 
noble liquid time projection chambers (TPCs) offer par- 
ticularly attractive prospects for dark matter detection 
through excellent background rejection capability, 3D po- 
sition sensitivity allowing definition of a fiducial volume, 
and cost-effective scalability [3"HTT]. Liquid xenon (LXe) 
in particular offers perhaps the most promising prospect 
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for unambiguous detection: it is intrinsically radio-pure, 
has the capability of powerful self-shielding p^HT4] . and 
is sensitive to low energy nuclear recoils. The background 
rejection power of two-phase TPCs comes from the ability 
to record both direct scintillation light (SI) and electro- 
luminescence from ionization (S2) of the target following 
an energy deposition. The ratio of the signal strength 
in these channels differs for electron and nuclear recoil 
interactions, allowing efficient discrimination between in- 
cident particle species. 

Traditionally the nuclear recoil energy scale for LXe 
TPCs with photomultiplier (PMT) read-out is set with 
the prompt primary scintillation channel SI. However, 
this channel is dominated by statistical fluctuations in 
signal generation and detection that ultimately preclude 
a direct mapping from observed SI to recoil energy. Al- 
though sufficient for setting limits on WIMP-nucleon 
cross-sections where no signal is observed, the SI energy 
scale severely limits the ability of a LXe TPC to test 
any possible observed signal excess against the WIMP 
hypothesis through analysis of the energy spectrum and 
event distribution. Furthermore, the requirement of SI 
signal across multiple PMTs imposes a relatively high 
energy threshold. In attempts to circumvent such limita- 
tions, an alternative approach is based on the ionization 
S2 channel alone. Although allowing a considerably lower 
threshold with sensitivity down to single electrons and 
a superior energy resolution due to electroluminescence 
gain in the gas phase, the ability to discriminate back- 
ground is lost and, in the absence of a prompt SI, the 
3D vertex identification is compromised. Furthermore, 
significant background (from, for e.g., VUV photoioniza- 
tion [TSJ [TB]) may limit the sensitivity of such S2-only 
analyses. 

In an effort to retain the discrimination power and fidu- 
cialization capability of the LXe TPC whilst simultane- 
ously exploiting the low energy threshold and superior 
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resolution afforded to the ionization channel, the adop- 
tion of a new data analysis method is proposed here, 
ft is shown that the traditional representation of data 
in Log(S2/Sl) vs. SI phase space can be replaced with 
Log(Sl/C2) vs. C2 (where C2 is the ionization yield) 
with no increase in background but with significant im- 
provement in energy resolution allowing reduction in the 
uncertainty in conversion to true recoil energy. This new 
technique, whilst retaining all the benefits of the LXe 
TPC technology such as self-shielding through fiducial 
definition and particle species discrimination, allows the 
direct conversion to an energy on an event-by-event basis. 
Proposed combined (scintillation and ionization) nuclear 
recoil energy scales [L7HT9] are still affected by the poor 
resolution in the SI channel for low energy recoils, whilst 
C2 represents a more appropriate energy scale. 

This article is organized as follows. In Section [TT] the 
traditional energy estimator, SI, and the ionization chan- 
nel, C2, are contrasted and modeled. Limitations of SI 
as the energy scale are highlighted, and C2 shown to pro- 
vide superior energy resolution and lower detector energy 
threshold for Xe. In Section 1111] the Log(S2/Sl) vs. SI 



and Log(Sl/C2) vs. C2 phase spaces are examined with 
the aid of a simple yet realistic model and monoenergetic 
neutrons. This illustrates the inefficiency in signal accep- 
tance, especially at higher energies, and misrepresenta- 
tion of energy in the traditional phase space. Appropri- 
ate signal acceptance regions, with clear energy thresh- 
olds that need not exploit smearing of sub-threshold en- 
ergy depositions into signal regions, are defined in the 
Log(Sl/C2) vs. C2 phase space. The Log(El/E2) vs. 
E2 phase-space is then introduced - these quantities are 
no longer detector specific and accessible as recoil energy 
for direct comparisons with other technologies. Finally, 
the response of the model detector to WIMPs is exam- 
ined in light of the new phase-space. 



II. NUCLEAR RECOIL ENERGY SCALE 



A. SI and C2 as nuclear recoil energy estimators 



liquid surface and travel through the gas phase, causing 
electroluminescence and giving a secondary scintillation 
signal (S2). The conversion from C2 to S2 is detector 
dependent - related to the gas gap and electric field con- 
figuration of the detector - it is simply the number of 
photoelectrons per electron extracted into the gas phase. 
For this reason C2, as the more meaningful parameter, is 
used as the primary observable rather than S2 through- 
out this article. 
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FIG. 1: Diagram of the production of excitation (SI) and ion- 
ization (C2, and subsequent S2) signals in 2-phase liquid/gas 
detectors (adapted from [20] ). 



Using SI as the primary nuclear recoil energy estimator 
appears to suffer from three fundamental difficulties: 

1. The SI signal is extremely small: typically a frac- 
tion of single photoelectron (PE) per keV at ener- 
gies below 10 keV nr . 

2. Conversion from SI to nuclear recoil energy, E nr 
depends on several parameters, given in Eq. [T] 



Fig. [T] illustrates the mechanism by which two-phase 
LXe TPCs record both the so-called SI and S2 signals 
following an interaction in the target volume. Energy de- 
position in LXe leaves atoms in both excited and ionized 
states. Excited Xe atoms combine with un-excited Xe 
atoms in a matter of a few picoseconds to form an excited 
Xe*, molecule. The de-excitation of this molecule releases 
ultraviolet photons. Simultaneously, ionized Xe atoms 
combine with non-ionized atoms to form XeJ molecules. 
In the absence of an electric field, these recombine with 
electrons to form an excited Xe** atom that emits ultra- 
violet photons through de-excitation (SI). In the pres- 
ence of an electric field, free electrons are not able to re- 
combine with the ionized Xe molecule and are drifted up 
to the liquid/gas interface. Under a strong enough elec- 
tric field, these free electrons (C2) are extracted from the 



SI 1 S n 



(i) 



where S m and S ee describe the electric field sup- 
pression factors for nuclear and electron recoils, re- 
spectively, C y is the detector light yield, and £ e ff 
is the scintillation yield for nuclear recoils. 

First the absolute C y calibration point is given by 
122 keV ee (keV electron equivalent) 57 Co 7-rays. 
Second, the difference between field dependent de- 
tector response to 7-rays and nuclear recoils must 
be accounted for (S nr and 5 CC ). Finally, the energy 
dependence of the scintillation yield for nuclear re- 
coils is corrected by £ e ff (shown in Fig. [2]) . 
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3. £ e ff is suppressed at lower energies, resulting in 
poorer energy measurements due to a lack of pho- 
toelectron statistics. 

It is well known that such extrapolation (from 122 
keV ee 7-rays to a few keV nr neutrons) and conversions 
tend to propagate large systematic errors at each stage. 
In other words, this is a so-called 'piece-by-piece' cali- 
bration, in opposition to ideal 'end-to-end' calibrations 
that tend to cancel out systematic errors. The investi- 
gation presented here indicates that the ionization yield, 
C2, can act as a much improved energy estimator for the 
following three reasons: 

1. The yield of C2 is an order of magnitude larger: 
typically ~5 electrons per keV nr . 

2. Conversion from C2 to nuclear recoil energy is 
straightforward: 



C2 
2/ 



(2) 



depending only on the one parameter Q y , the 
charge yield, such that the systematic error in the 
conversion is reduced. The energy dependence of 
Q y is shown in Fig. [3j 

3. Theoretical models and dedicated measurements 
indicate that the value of Q y increases with de- 
creasing energy, meaning that even at low energy, 
large C2 signals are expected. 
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FIG. 2: Collection of all direct measurements of C e g to 
date [211425] . In this diagram, the solid grey regions repre- 
sent a Maximum Likelihood fit to all direct measurements 
and the associated la and 2a confidence bands; this parame- 
terization is adopted in this work and is taken from [SJ. £ e g 
is extrapolated to zero at 1 keV nr - In addition, a contempo- 
raneous indirect (data vs. Monte Carlo) measurement [26] is 
shown to be in good agreement with direct measurements. 
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FIG. 3: Charge yield (Qy) as measured directly at several 
electric fields |20l I27] . Also included are theoretical Q y curves 
derived using stopping powers as calculated by Hitachi [28] 
(dashed curve) and Lindhard [29] (solid curve) and extracted 
from [30]. For the purposes of this paper, a Q y that falls 
between the Hitachi and Lindhard measurements is adopted 
(red curve). 



B. Modeling SI and C2 signals 

In order to perform a further systematic comparison, 
a model is constructed using Monte Carlo simulations. 
The underlying, well-motivated assumptions [231 l2"ol |2"51 - 
I3"T] are as follows: 

• The energy dependence of C e s and Q y is as shown 
in Figs.[2]&[3j respectively. 

• £ y = 3 PE/keV C c at 122 keV C0 . 

• S nl and S^e are chosen assuming an electric field of 
0.5 kV/cm. 

• For conversion from S2 to C2, the secondary scin- 
tillation gain in a gas phase is assumed to be 20 
(i.e., S2 = C2 x 20). 

• Poisson distribution for SI in terms of the number 
of photoelectrons. 

• Poisson distribution for C2 in terms of the number 
of free electrons from initial ionization. 

At this stage new terminology may be introduced: 'El' 
as the energy estimated by SI, and 'E2' as the energy 
estimated by C2. The dominant source of the energy 
mis-measurement is the Poisson fluctuations in SI in the 
case of El, and in C2 for the case of E2. It should be 
noted that any other factors such as the secondary scin- 
tillation process from C2 to S2 have a negligible contri- 
bution. E2 by itself represents a better estimate of the 
energy than the combination of El and E2 (as adopted 
for electromagnetic interactions of »10 keV ee ) due to the 
poor photoelectron statistics of SI. 

Fig. [4] shows the simple relationship between SI, C2 
and nuclear recoil energy. Here three different cases, 
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C y = 2, 3, & 5 PE/keV C o, are shown for SI. For a given 
energy, C2 is an order of magnitude larger than SI, and 
the advantage of C2 over SI grows for lower energies. 
Fig. [5] shows the energy resolution (c/E) of El and E2. 
Here again the three cases of C y — 2, 3, & 5 PE/keV co 
are presented. Even with C y as large as 5 PE/keV cc , El 
becomes severely limited below ~10 keV nr . On the other 
hand, E2 shows reasonable resolution (~20%) down to 3 
keVnr- 
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FIG. 4: SI, for three different detector light yields, and C2 
observables as a function of nuclear recoil energy. SI is shown 
in units of PE, whereas C2 is in electrons (i.e., prior to am- 
plification to S2). 




FIG. 5: Energy resolution as a function of recoil energy for 
an energy scale determined using SI (for three different light 
yields) and C2. Above ~40 keV nr , the resolutions become 
comparable but at lower recoil energies, of relevance to WIMP 
searches, the advantage of a C2 derived energy scale is evident. 



III. COMPARISON OF THE DIFFERENT 
PHASE SPACES 

To visualize the problem of the conventional method, 
Log(S2/Sl) vs. SI, four phase spaces are contrasted: 

I. S2 vs. SI 

II. Log(S2/Sl) vs. SI 

III. Log(Sl/C2) vs. C2 

IV. Log(El/E2) vs. E2 

Recall that SI is analytically connected to El by £ e ff, 
while C2 is analytically connected to E2 by Q y . Thus 
there exists a strict one-to-one mapping between all four 
phase spaces listed. 



A. Monoenergetic nuclear recoil injection 

As a simple illustration of these phase spaces, we first 
inject monoenergetic nuclear recoil energies of 4, 8, 16 
and 32 keV nr . Such energies are representative of the 
range of interest for WIMP scatters in LXe. Typically de- 
tectors are calibrated with broad spectrum sources such 
as Am-Be or 252 Cf neutron emitters that produce recoils 
across this energy range. Note that the results are of 
course unaffected for the case of continuous spectra. 

Firstly, Fig. [6] shows S2 vs. SI. Here, as well as 
throughout the article, the following are indicated: 

• For each of the sets of monoenergetic energy recoils, 
lcr and 2a regions are shown as contours. 

• A black solid curve indicates E2 = El. Additional 
±lcr and ±2cr curves show the band of continuous 
nuclear recoil energy distribution (as from, for e.g., 
Am-Be neutron calibration data). 

• El = 4, 8, 16, 32 keV nr are indicated with green 
lines, while E2 = 4, 8, 16, 32 keV nr are indicated 
with blue lines. 

It is clearly seen that at lower energy the SI distribu- 
tion is much wider than the S2 distribution, highlight- 
ing the fact that S2 is the better energy estimator. At 
4 keV nr , SI is only 2 PE, thus its Poisson fluctuation 
prohibits any useful energy estimate. 

Next, Fig. [7] shows the traditional phase space of 
Log(S2/Sl) vs. SI for the monoenergetic recoils. It is 
seen that the recoils are distributed in an elongated el- 
lipse for any given energy. The ellipses are aligned with 
contours of constant E2 because E2 represents the nu- 
clear recoil energy with less spread (and ambiguity) than 
SI. It is evident that this phase space is distorted, and 
thus the X-axis does not represent the energy of events. 

Fig. 8 introduces the new phase space of Log(Sl/C2) 
vs. C2. Unlike in the Log(S2/Sl) vs. SI phase space 
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(Fig. [7]) , each energy band is now well aligned as a ver- 
tical ellipse. Note that they are still well separated, 
as in Fig. [6] showing S2 vs. SI. Note also the sim- 
ple relation Log(Sl/C2) = Log(Sl/(S2/20)) = 1.30 - 
Log(S2/Sl). Therefore the Y-axis of this new phase 
space, Log(Sl/C2), is merely a mathematical transfor- 
mation of Log(S2/Sl). The critical advantage of this 
method comes from adopting the new X-axis of C2 in- 
stead SI, whilst retaining discrimination. 



E1^4keV nr SkeV m 16keV M 32 keV n 




FIG. 6: The relationship between SI and S2 for monoener- 
getic nuclear recoils of 4, 8, 16, and 32 keV nr . Solid contours 
represent la (blue) and 2a (grey) deviations from the mean 
for these monoenergetic neutrons. Dashed curves represent 
the la (blue) and 2a (grey) deviations from the mean from 
continuous energy recoils. Green lines (vertical here) indicate 
energy in El and blue lines (horizontal here) indicate energy 
in E2. Note that the same convention is used throughout this 
article. It is clear that the Poisson fluctuations in SI gives a 
greater distortion for lower recoil energies making an energy 
scale determined by SI significantly more ambiguous as true 
recoil energy gets lower. 



Finally, in Fig. [9j we introduce the phase space of 
Log(El/E2) vs. E2: 

• The X-axis directly represents the best estimate of 
the nuclear recoil energy on an event-by-event ba- 
sis. Fig. [P0| further highlights the superior energy 
resolution of the E2 phase space. 

• Log(El/E2) = (i.e., E2 = El) sets the mean for 
all the nuclear recoil events, in the absence of im- 
posed thresholds. 

• Experimental specific calibrations are incorporated 
into the only physically meaningful parameters, E2 
(energy) and E1/E2 (discrimination); this phase 
space is not specific to any single detector. 



El =4keV„8keV n , ldkeV n , 
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FIG. 7: The distribution of monoenergetic nuclear recoils as 
shown in traditional Log(S2/Sl) vs. SI parameter space. The 
distortion of what has been previously used to represent a 
recoil energy scale is seen. The SI is not a well constrained 
representative of energy on an event-by-event basis. 
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FIG. 8: The monoenergetic nuclear recoils represented in the 
Log(Sl/C2) vs. C2 phase space. Transposed onto the C2 
scale, the ambiguity of the conversion between an observable 
and a recoil energy is significantly reduced, with no loss of 
discrimination. 



B. Representing the traditional signal box 

To further evaluate the benefits and disadvantages of 
the four different phase spaces, let us now consider a 
'signal box' for a WIMP search. Traditionally the WIMP 
signal box is defined by: 

• The accepted energy range, determined by a mini- 
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E2 = 4 keV^ 8 keV„ ISkeV^ 32 keV„ 
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FIG. 9: For the comparison of observed and expected WIMP 
recoil signals, it is necessary to translate into a true nuclear 
recoil energy scale, as represented here in Log(El/E2) vs. E2 
phase space. 



mum and maximum energy of El (given by SI). 

• The neutron band given by Am-Be calibration 
data. In Log(S2/Sl), events below the mean of 
the distribution (where discrimination is highly ef- 
ficient) and within la or 3a are accepted. 

As a typical and realistic example, let us assume: 

• An El range of 8 - 32 keV nr 




E2 [keVnr] 



• The neutron population in the band between the 
mean and 2a in Log(S2/Sl) vs. SI. 



This traditional signal box is shown in Fig. 11 as con- 
tained within the pink borders in all four phase spaces. 
There are two important remakes to be made: 



1 



The maximum energy of 32 keV nr is indicated in 
El as a vertical line in the traditional Log(S2/Sl) 
vs. SI phase space. However, only a very small 
fraction of 32 keV nr events are actually accepted 



by the signal box (see Fig. 13). 



2. At the minimum energy of 8 keV nr , the bound is 
rather loose towards low energy events. There are 
some fraction of events between 4 and 8 keV nr that 
penetrate into the signal box. This is the origin of 
the so-called 'Poisson smearing' that enhances the 
sensitivity of such detectors to low energy events, 
and thus to low mass (<10 GeV/c 2 ) WIMPs whilst 
maintaining a relatively high threshold (8 keV nr in 
this model) . Here some fraction of events that de- 
posit energies below threshold may be detected in 
the signal region following statistical fluctuations 
in the number of detected PE. 



FIG. 10: Comparison of 8 keV nr and 16 keVnr nuclear recoil 
energy spectra in El (top) and E2 (bottom), highlighting the 
superior energy resolution in E2 that improves relative to El 
for falling energies. 



It is worth noting that the traditional signal region 
shows up as quite a distorted box in both Log(Sl/C2) 



vs. C2 and Log(El/E2) vs. E2, as shown in Fig. 11 



C. Introduction of a new signal box 

To address the fact that El is not a good energy esti- 
mator, a new signal box, based on E2, is defined. In the 
case of E2, the energy can be well determined even down 
to the 2-3 keV nr range. Therefore, there is no reason to 
define a lower bound at 8 keV nr (which is necessary in the 
case of the El-based phase space to ensure good signal 
acceptance at threshold). Instead, a minimum energy of 
4 keV nr (close to the lowest directly measured relative 
scintillation and ionization yield values) is introduced in 
this model. The signal box that is proposed is shown in 
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FIG. 11: The traditional signal box (pink boundaries) as denned in Log(S2/Sl) vs. SI and applied to the four phase spaces 
described in the text. It can be seen that WIMP recoil acceptance at 4 and 32 keV nr is severely reduced in relation to that 
suggested in the SI space once the box is transposed into the C2 or E2 energy scales. 



Fig. [12] as a red box in all four phase spaces and is defined 
by: 

• E2 range between 4 and 32 keV nr - 

• The region between the mean and 2er of the neutron 
band in Log(Sl/C2) vs. C2 or Log(El/E2) vs. E2. 

To further illustrate the effect of minimum and max- 
imum energy bounds in the old and new methods, the 
signal acceptance is shown as a function of true recoil 
energy in Fig. [13] The new signal box shows a flat ac- 
ceptance in energy from 4-32 keV nr , whilst an El-based 
signal region compromises efficiency both at low and high 
energies. 

As a final extension to the issue of signal box definition, 
it should be noted that in the phase space of Log(El/E2) 
vs. E2, a sharp lower energy bound at 4 keV nr is unnec- 
essary. Assuming extrapolation of models for relative 
scintillation and ionization yields to lower energies, or in 



light of future direct measurements elucidating this re- 
gion, a lower bound that follows the minimum contour 
of 2 PE in the SI channel may be adopted. Similarly, 
there is no reason to set a higher energy bound. Indeed, 
a profile likelihood analysis across the 2D phase space 
of Log(El/E2) vs. E2 allows improved background re- 
jection and signal excess energy distribution analysis as 
a result of the improved resolution. Such a maximally 
efficient signal region is presented in Fig. |14| 



D. WIMP spectra injection 

As the final stage of this investigation, injection of 
WIMP energy spectra are simulated. For illustra- 
tion purposes a WIMP mass of 100 GeV/c 2 is as- 
sumed with a spin-independent WIMP-nucleon interac- 
tion cross-section of 5xl0~ 45 cm 2 . An exposure of 1 year 



for a 1 ton fiducial mass is considered. Fig. 15 shows 
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FIG. 12: A new signal box (red boundaries) is defined in Log(El/E2) v E2 phase space and applied to the four phase spaces 
described in the text. Applying this box into traditional phase space of Log(S2/Sl) vs. SI serves to highlight the huge distortion 
present in the traditional method of signal box definition and analysis. The C2 and E2 energy scales also allow measurements 
of discrimination parameters and estimation of background leakage into the signal region from slices in C2 specific to a given 
energy. 



the WIMP recoil event distributions in the traditional 
Log(S2/Sl) vs. SI phase space and signal box as pre- 
viously defined in the SI energy estimator scale (as in 
Fig. 111. Fig. 16 shows the same WIMP event distri- 
bution in Log(El/E2) vs. E2 with the signal box as 
defined with a lower threshold dictated only by the 2 
PE requirement for SI (as in Fig. 



14). The lower en- 



ergy threshold from the E2 scale and signal box bounded 
only by the SI threshold naturally results in significantly 
greater signal acceptance for WIMPs: of order 500 events 
are observed as compared to approximately 300 for the 
traditional box. 

Furthermore, with the E2 energy scale, given appropri- 
ate statistics, the new representation can expose an expo- 
nential energy spectrum with a WIMP mass-dependent 
slope with greater confidence than the traditional SI en- 
ergy scale. This has an important consequence on the 



discovery of WIMPs. The new method of Log(El /E2) vs. 
E2 allows investigation of not only the excess of events 
over expected background, but also the energy spectrum 
of these events and determination of WIMP mass. 



IV. DISCUSSION AND CONCLUSIONS 
A. Discussion 

The new concept of Log(El/E2) vs. E2 as a replace- 
ment of the current standard of Log(S2/Sl) vs. SI has 
been introduced based on Monte Carlo simulations of 
a theoretically motivated and detector realistic model. 
However, the true power of the new analysis method can 
be realized only in conjunction with real data analysis. 
In particular, the method makes for a considerably more 
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An Ideal Sig. Box (S 1 thr. = 2 PE only) 



E2-based Sig. Box (E2 = 4 - 32 keVnr) 

Traditional El-based Sig. Box (El - 8 - 32 keVnr) 
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FIG. 13: Signal acceptance in case of the traditional and 
newly proposed signal boxes. The El (SI) based acceptance 
(solid purple curve) when defined as 8 - 32 keV nr in SI is in- 
efficient at higher energies when compared to the same upper 
bound as defined in E2 (C2) (dashed red curve). The lower 
threshold afforded this signal box increases the efficiency at 
lower energies. Indeed, contingent upon the behavior of elec- 
tromagnetic background, in principle the lower bound may be 
completely removed and the threshold defined simply as the 
minimum requirement of 2 PE (i.e., 1 PE in 2 PMTs). No 
upper bound need be enforced at all. This is represented by 
the solid red curve. 



stringent probe of accurate reproduction of neutron cal- 
ibration data by Monte Carlo simulations. Demanding 
coherence across all four types of phase space, systematic 
comparison between data and Monte Carlo simulations 
will allow considerably more robust estimates of nuclear 
recoil acceptance, as well as verification of assumptions 
in the energy dependence of £ G ff and Q y . 

One should note, however, that the most reliable tests 
of such models could be performed by direct calibration 
of dark matter detectors in situ, using monoenergetic 
neutron generators, very much in the same manner as 
dedicated £ e ff measurements. 

With well modeled neutron calibration data, repro- 
duced by simulation across all phase spaces, the detection 
efficiency can be reliably measured down to very low en- 
ergy. Provided background remains negligible, the low- 
est possible energy threshold in both SI and C2 can be 
adopted. Although for C2 this could be as low as several 
electrons (~keV nr ), assuming sufficient electron mobility 
in the liquid, the minimum requirement of SI = 2 PE 
nonetheless dominates. 

The proposed approach also has consequence for back- 
ground expectation analysis. In reality, since the energy 
may be determined more reliably with the new method, 
the background distribution in the new phase space (ei- 
ther Log(Sl/C2) vs. C2 or Log(El/E2) vs. E2) is 
more meaningful. The traditional approach of estimating 



FIG. 14: The upper and lower energy bounds in E2 that 
are present in Fig. [12] are shown as red dashed vertical lines. 
These may be relaxed such that the signal region is defined 
by the solid red boundaries. Here no upper boundary is im- 
posed, and the lower boundary is defined by the SI = 2 PE 
coincidence requirement. Acceptance to WIMP recoils is thus 
enhanced. This signal acceptance for this box is the corre- 
sponding solid red curve of Fig. [13] 




FIG. 15: WIMP distribution in traditional phase space of 
Log(S2/Sl) vs. SI (assuming 100 GeV/c 2 WIMP mass, 
WIMP-nucleon interaction cross-section of 5xl0~ 45 cm 2 ) for 
a 1 ton fiducial volume and 1 year exposure. A total of 277 
events are observed in the defined signal region. 



background in slices of fixed SI is misleading because in 
each such slice, the energy of events is not constant; the 
larger the value of (S2/S1), the larger the energy. Fur- 
thermore, the new phase space may be of considerable 
benefit in characterising rare topology background, such 
as so-called Multi-Scintillation Single-Ionization (MSSI) 
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FIG. 16: WIMP distribution in the proposed phase space 
of Log(El/E2) vs. E2 (assuming 100 GeV/c 2 WIMP mass, 
WIMP-nucleon interaction cross-section of 5xl0~ 45 cm 2 ) for 
a 1 ton fiducial volume and 1 year exposure. A total of 513 
events are observed in the defined signal region. 



events. These are predominantly 7-ray multiple scatter 
events where one vertex is in the bulk volume, produc- 
ing both scintillation (SI) and ionization (S2), but the 
second is in a region with no applied drift field (such as 
below the TPC cathode). Partial scintillation (SI', El') 
may be detected from such 'inactive' volumes whereas 
no electrons are extracted into the gas phase meaning 
no ionization signal. Such events appear with artificially 
large SI signals (from both interaction sites) relative to 
the ionization signal (from a single site) and leak towards 
the neutron band. In the traditional phase space this 
causes a distortion in both the Y-axis, Log(S2/Sl), and 
X-axis, SI, such that: 

[Log(S2/Sl), Sl]-[Log(S2/(Sl + Sl')), (Sl + Sl')] 

In the new space only the Y-axis is affected: 

[Log(El/E2), E2] -» [Log((El + El')/E2), E2] 

The result is that the MSSI events produce a flat distri- 
bution in the E2 phase space that is more readily char- 
acterised. 

The consequences of the new approach on the elec- 
tron recoil rejection is being explored by developing a 
similar model of (SI, C2) observable signals from 7-rays. 
However it is clear that such an approach (translating to 
different phase spaces) cannot introduce any additional 
background events to a signal box - every event in any of 
the four phase spaces can be mapped to any other phase 
space. The number of background events in one phase 
space remains exactly the same in any other. 

Finally, with regards any future positive signal, since 
the proposed method dramatically improves both the 



acceptance and the energy determination, the WIMP 
mass can be determined far more reliably than as pre- 
sented in previous studies for the future large-scale de- 
tectors [32 , 33J . A thorough exploration with this revised 
analysis will be presented in a separate article. 



B. Conclusions 

Traditionally, nuclear recoil energy for WIMP searches 
by liquid Xe dual-phase TPCs has been estimated by 
means of the initial scintillation and not by the ionization 
charge. The secondary scintillation is exploited only as 
a discrimination parameter in Log(S2/Sl) vs. SI. Alter- 
native approaches that use only the S2 channel seriously 
damage the discrimination power and 3D event vertex 
information. 

This article has presented issues with the traditional 
approach, primarily due to the fluctuations in SI signal 
that do not allow reliable determination of the true re- 
coil energy. Additionally, the acceptance when defined 
in SI does not translate to the assumed nuclear recoil 
energy, causing efficiency loss even towards the higher 
energy boundary of the signal box. 

To address these problems, the new phase space of 
Log(Sl/C2) vs. C2 has been introduced using C2 as 
the superior energy estimator. Crucially, the improved 
resolution can be exploited without loss of fiducializa- 
tion (critical to liquid noble gas TPCs for self-shielding 
and event reconstruction) or discrimination (critical in 
rejecting single scatter electronic background). 

Finally, two energy estimators are defined (El from 
SI, and E2 from C2), and Log(El/E2) vs. E2 is intro- 
duced. Advantages of this approach may be summarized 
as follows: 

1. The energy resolution from E2 is considerably en- 
hanced relative to the El, especially at low energies 
where the WIMP spectrum peaks. This provides 
improvement in evaluating the WIMP energy spec- 
trum against the background. 

2. Lower energy thresholds and greater signal accep- 
tance can be achieved without compromizing self- 
shielding or discrimination power. 

3. Monte Carlo simulation models of neutron calibra- 
tions can be assessed in multiple phase spaces to 
address underlying assumptions and nuclear recoil 
efficiency can be uniquely extracted. 

4. For a given E2 bin, Log(El/E2) represents an elec- 
tronic rejection parameter specific to that energy. 

5. With multiple WIMP events observed the WIMP 
mass can be estimated with greater confidence. 

In conclusion, the proposed method of Log(El/E2) vs. 
E2 is more accurate and reliable than the traditional 
method of Log(S2/Sl) vs. SI. El and E2 are presented 
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as more accessible detector independent parameters that 
can be adopted within the dark matter search community 
and directly compared with other technologies. 
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